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^ This is a summary of presentations delivered at the OCl parallel session "Primordial 

, Gravitational Waves and the CMB" of the 12**^ Marcel Grossmann meeting in Paris, July 

' 2009. The reports and discussions demonstrated significant progress that was achieved in 

theory and observations. It appears that the existing data provide some indications of the 
presence of gravitational wrave contribution to the CMB anisotropies, while ongoing and 
planned observational efforts are likely to convert these indications into more confident 
statements about the actual detection. 
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1. Introduction 

^ , Here, we report on presentations and discussions that took place at the OCl parallel 

I session entitled "Primordial Gravitational Waves and the CMB" . The programme 

OO ' was designed to include both theoretical and observational results. The electronic 

I versions of some of the talks delivered at this session are available at the MG12 

' website.'i' We shall start from the overall framework of this session and motivations 

for studying the primordial (relic) gravitational waves. 

Primordial gravitational waves (PGWs) are necessarily generated by a strong 
variable gravitational field of the very early Universe.I^ The existence of relic gr- 
tavitational waves relies only on the validity of basic laws of general relativity and 
quantum mechanics. Specifically, the generating mechanism is the superadiabatic 
■ (parametric) amplification of the waves' zero-point quantum oscillations. In con- 

trast to other known masslcss particles, the coupling of gravitational waves to the 
external ( "pump" ) gravitational field is such that they could be classically ampli- 
fied or quantum- mechanically generated by the gravitational field of a homogeneous 
isotropic FLRW (Friedmann-Lemaitre-Robertson- Walker) universe. Under certain 
extra conditions the same applies to the primordial density perturbations. The 
PGWs are the cleanest probe of the physical conditions in the early Universe right 
down to the limits of applicability of currently available theories, i.e. the Planck den- 
sity ppi = c'/G'^n « lO^-^g/cm^ and the Planck size h\ = {Gh/c^y/'^ « lO'^^cm. 

The amount and spectral content of the PGWs field depend on the evolution 
of the cosmological scale factor a{T]) representing the gravitational pump field. The 
theory was applied to a variety of a{ri), including those that are now called infla- 
tionary models.'^^ If the exact a{r]) were known in advance from some fundamental 
"theory-of-everything" , we would have derived the properties of the today's signal 
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with no ambiguity. In the absence of such a theory, we have to use the available 
partial information in order to reduce the number of options. The prize is very high 
- the actual detection of a particular background of PGWs will provide us with a 
unique clue to the birth of the Universe and its very early dynamical behaviour. 

To be more specific, let us put PGWs in the context of a complete cosmological 
theory hypothesizing that the observed Universe has come to the existence with 
near-Planckian energy density and size (see paper and references therein). It 
seems reasonable to conjecture that the embryo Universe was created by a quantum- 
gravity or by a "theory-of-everything" process in a near-Planckian state and then 
started to expand. (If you think that the development of a big universe from a tiny 
embryo is arrant nonsense, you should recollect that you have also developed from a 
single cell of microscopic size. Analogy proposed by the biophysicist E. Grishchuk.) 
The total energy, including gravity, of the emerging classical configuration was likely 
to be zero then and remains zero now. 

In order for the natural hypothesis of spontaneous birth of the observed Universe 
to bring us anywhere near our present state characterized by pp = SHq/SttG ~ 
10~^^g/cm^ and Ip = 5001h (which is the minimum size of the present-day patch of 
homogeneity and isotropy, as follows from observation^^ the newly-born Universe 
needs a significant 'primordial kick'. During the kick, the size of the Universe (or, 
better to say, the size of our patch of homogeneity and isotropy) should increase 
by about 33 orders of magnitude without losing too much of the energy density of 
whatever substance that was there, or maybe even slightly increasing this energy 
density at the expence of the energy density of the gravitational field. This process 
is graphically depicted in Fig. [1] (adopted from the pa.pe^. The present state of the 
accessible Universe is marked by the point P, the birth of the Universe is marked by 
the point B. If the configuration starts at the point B and then expands according to 
the usual laws of radiation-dominated and matter-dominated evolution (blue curve) , 
it completely misses the desired point P. By the time the Universe has reached the 
size Ip, the energy density of its matter content would have dropped to the level 
many orders of magnitude lower than the required pp. The only way to reach P 
from B is to assume that the newly-born Universe has experienced a primordial 
kick allowing the point of evolution to jump over from the blue curve to the black 
curve. 

If we were interested only in the zero-order approximation of homogeneity and 
isotropy, there would be many evolutionary paths equally good for connecting the 
points B and P. However, in the next-order approximations, which take into account 
the inevitable quantum-mechanical generation of cosmological perturbations, the 
positioning and form of the transition curve in Fig.[T]become crucial. The numerical 
value of the Hubble parameter H (related to the energy density of matter driving 
the kick, as shown on the vertical axis of the figure) determines the numerical level 
of amplitudes of the generated cosmological perturbations, while the shape of the 
transition curve determines the shape of the primordial power spectrum. 

The simplest assumption about the initial kick is that its entire duration was 
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Fig. 1. A primordial kick is required in order to reach the present state of the Universe P from the 
birth event B. Red lines describe possible transitions that would be accompanied by the generated 
cosmological perturbations of observationally necessary level and spectral shape0The "legitimate 
transition in inflationary theory" is an evolution allowed by the incorrect (inflationary) formula 
for density perturbations. See explanations in Sec l3.2l below. 

characterized by a single power-law scale factoil^l 

a(77)-?ohr+^ (1) 

where Iq and /? are constants, and /3 < — 1. In this power-law case, the gravitational 
pump field is such that the generated primordial metric power spectra (primordial 
means considered for wavelengths longer than the Hubble radius at the given mo- 
ment of time), for both gravitational waves and density perturbations, have the 
universal power-law dependence on the wavenumber n: 

h\n) oc (2) 

It is common to write these metric power spectra separately for gravitational waves 
(gw) and density perturbations (dp): 

h'ln) (gw) = B^n"', h'in) (dp) = Bfn"^-!. (3) 

According to the theory of quantum-mechanical generation of cosmological 
perturbations (for a review, se^, the spectral indices are approximately equal, 
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Us — 1 = rit = 2{/3 + 2), and the amplitudes Bt and Bg are of the order of magni- 
tude of the ratio Hi/Hp\, where Hi ^ c/lo is the characteristic value of the Hubble 
parameter H during the kick. The straight lines in Fig. [T] symbolize the power-law 
kicks (HI) . They generate primordial spectra with constant spectral indices through- 
out all wavelenghts. In particular, any horizontal line describes an interval of de 
Sitter evolution, f3 ~ —2, H — 0, H — const. (Initial kick driven by a scalar field 
is appropriately called inflation: dramatic increase in size with no real change in 
purchasing power, i.e. in matter energy density.) The gravitational pump field of 
a de Sitter kick transition generates perturbations with flat (scale-independent) 
spectra Us ~ 1 — nt = 0. The red horizontal line shown in Fig. [T] corresponds to 
Hi/Hpi « 10^^ and the generated primordial amplitudes Bt ~ Bs ~ 10~^. In nu- 
merical calculations, the primordial metric power spectra are usually parameterized 

by 

Pt{k) = At{k/kor\ P,(k)^As{k/koT^-\ (4) 

where fco = 0.002Mpc-^ 

Although the assumption of a single piece of power-law evolution ([T]) is simple 
and easy to analyze, the reality could be more complicated and probably was more 
complicated (see the discussion of WMAP data in Sec l3.2[) . A less simplistic kick can 
be approximated by a sequence of power-law evolutions, and then the primordial 
power spectra will consist of a sequence of power-law intervals. 

The amplitudes of generated cosmological perturbations are large at long wave- 
lengths. According to the widely accepted assumption the observed anisotropies in 
the cosmic microwave background radiation (CMB) are caused by perturbations 
of quantum-mechanical origin. This assumption is partially supported by the ob- 
served "peak and dip" structure of the CMB angular power spectra. Presumably, 
this structure reflects the phenomenon of quantum-mechanical phase squeezing and 
standing- wave pattern of the generated metric fields.'^ The search for the relic grav- 
itational waves is a goal of a number of current and future space-borne, sub-orbital 
and ground-based CMB experiments.!^^ 

The CMB anisotropies are usually characterized by the four angular power spec- 
tra Cj'^ ^ Cf^, Cf^ and Cj^ as functions of the multipole number i. The contri- 
bution of gravitational waves to these power spectra has been studied, both analyt- 
ically and numerically, in a number of papers.^i^ll^II xhe derivation of today's CMB 
power spectra brings us to approximate formulas of the following structurePSl 

Cr = r-iHn,lrec)? [Fl{n)\\ 

Cr = I ^Hn, Vrec)h'{n, r^rec) [Fl{n)Ff{n)] , (5) 
Cr = I ^[h'{n,T^rec)? [Fl{n)]\ where Y^E^B. 

In the above expressions, [h{n,rirec)]'^ and [h'{n,r]rec)]^ are power spectra of the 
gravitational wave field and its first time-derivative. The spectra are taken at the 
recombination (decoupling) time r/rec- The functions {n) {X = T, E, B) take care 
of the radiative transfer of CMB photons in the presence of metric perturbations. To 
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a good approximation, the power residing in the metric fluctuations at wavenumber 
n translates into the CMB TT power at the angular scales characterized by the 
multipolcs £ ^ n. Similar results hold for the CMB power spectra induced by 
density perturbations. Actual numerical calculations use equations more accurate 
than Eq. dS]). 

There are several differences between the CMB signals arising from density 
perturbations and gravitational waves. For example, gravitational waves produce 
B-component of polarization, while density perturbations do notf^^- gravitational 
waves produce negative TE-correlation at lower multipoles, while density perturba- 
tions produce positive correlation,'2l2ni2lH2l] g^j^^ gQ However, it is important to 
realize that it is not simply the difference between zero and non-zero or between 
positive and negative that matters. (In any case, various noises and systematics 
will make this division much less clear cut.) What really matters is that the gw 
and dp sources produce different CMB outcomes, and they are in principle distin- 
guishable, even if they both are non-zero and of the same sign. For example, if the 
parameters of density perturbations were precisely known from other observations, 
any observed deviation from the expected TT, TE and EE correlation functions 
could be attributed (in conditions of negligible noises) to gravitational waves. From 
this perspective, the identification of the PGWs signal goes well beyond the often 
stated goal of detecting the B-mode of polarization. In fact, as was argued in the 
talk by D. Baskaran, delivered on behalf of the group including also L. P. Grishchuk 
and W. Zhao, the TT, TE, and EE observational channels could be much more 
informative than the BB channel alone. Specifically for the Planck mission, the in- 
clusion of other correlation functions, in addition to BB, will significantly increase 
the expected signal-to-noise ratio in the PGWs detection. 

It is convenient to compare the gravitational wave signal in CMB with that 
induced by density perturbation. A useful measure, directly related to observations, 
is the quadrupole ratio R defined by 

i.e. the ratio of contributions of gw and dp to the CMB temperature quadrupole. 
Another measure is the so-called tensor-to-scalar ratio r. This quantity is built from 
primordial power spectra (0]): 

Usually, one finds this parameter linked to incorrect (infiationary) statements. Con- 
cretely, infiationary theory substitutes (for "consistency" ) its prediction of arbitrar- 
ily large amplitudes of density perturbations in the limit of models where spectral 
index Ug approaches Ug — 1 and nt approaches nt — Q (horizontal transition lines in 
Fig. [1]) by the claim that it is the amount of relic gravitational waves, expressed in 
terms of r, that should be arbitrarily small. (For more details, see Sec. 13.21 ) How- 
ever, if r is defined by Eq.([7]) without implying infiationary claims, one can use this 
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parameter. Moreover, one can derive a relation between R and r which depends on 
the background cosmological model and spectral indices. For a rough comparison 
of results one can use r w 2R. 

2. Overview of oral presentations 

The OCX session opened with the talk of Brian Keating, appropriately entitled "The 
Birth Pangs of the Big Bang in the Light of BICEP" . The speaker reported on the 
initial results from the Background Imaging of Cosmic Extragalactic Polarization 
(BICEP) experiment. The conclusions were based on data from two years of ob- 
servation. For the first time, some meaningful limits on r were set exclusively from 
the fact of (non)observation of the B-mode of polarization (see figure [2] adopted 
from Chiang et al. papei'^S). The author paid attention to various systematic effects 
and showed that they were smaller than the statistical errors. B. Keating explained 
how the BICEP 's design and observational strategies will serve as a guide for fu- 
ture experiments aimed at detecting primordial gravitational waves. In conclusion, 
it was stressed that a further 90% increase in the amount of analyzed BICEP data 
is expected in the near future. 

The next talk, delivered by Deepak Baskaran (co-authors: L. P. Grishchuk and 
W. Zhao), was entitled "Stable indications of relic gravitational waves in WMAP 
data and forecasts for the Planck mission" . D. Baskaran reported on the results of 
the likelihood analysis, performed by this group of authors, of the WMAP 5-year 
TT and TE data at lower multipoles. Obviously, in the center of their effort was 
the search for the presence of a gravitational wave signal.l22l2il Por the parameter R, 
the authors found the maximum likelihood value R — 0.23, indicating a significant 
amount of gravitational waves. Unfortunately, this determination is surrounded by 
large uncertainties due to remaining noises. This means that the hypothesis of no 
gravitational waves, R — 0, cannot be excluded yet with any significant confidence. 
The speaker compared these findings with the result of WMAP team, which found 
no evidence for gravitational waves. The reasons by which the gw signal can be 
overlooked in a data analysis were identified and discussed. Finally, D. Baskaran 
presented the forecasts for the Planck mission. It was shown that the stable in- 
dications of relic gravitational waves in the WMAP data are likely to become a 
certainty in the Planck experiment. Specifically, if PGWs are characterized by the 
maximum likelihood parameters founcP^HH from WMAP5 data, they will be de- 
tected by Planck at the signal-to-noise level S/N = 3.65, even under unfavorable 
conditions in terms of instrumental noises and foregrounds. (For more details along 
these lines, see sections below.) 

The theoretical aspects of generation of CMB polarization and temperature 
anisotropics by relic gravitational waves were reviewed by Alexander Polnarev in 
the contribution "CMB polarization generated by primordial gravitational waves. 
Analytical solutions" . The author described the analytical methods of solving the 
radiative transfer equations in the presence of gravitational waves. This problem 
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Fig. 2. BICEPs TE, EE, and BB power spectra together with data from other CMB polarization 
experiments. Theoretical spectra from a CDM model with r = 0.1 are shown for comparison. The 
BB curve is the sum of the gravitational wave and lensing components. At degree angular scales 
BICEPs constraints on BB are the most powerful to date.^SU 



is usually tackled by numerical codes, but this approach foreshadows the ability 
of the researcher to understand the physical origin of the results. The analytical 
methods are a useful complement to numerical techniques. They allow one not only 
to explain in terms of the underlying physics the existing features of the final prod- 
uct, but also anticipate the appearance of new features when the physical conditions 
change. A. Polnarev showed how the problem of CMB anisotropics induced by grav- 
itational waves can be reduced to a single integral equation. This equation can be 
further analyzed in terms of some integral and differential operators. Building on 
this technique, the author formulated analytical solutions as expansion series over 
gravitational wave frequency. He discussed the resonance generation of polarization 
and possible observational consequences of this effect. 

An overview of the current experimental situation was delivered by Carrie Mac- 
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Tavish in the talk "CMB from space and from a balloon". C. MacTavish focused 
on the interplay between experimental results from two CMB missions: the Planck 
satellite and the Spider balloon experiment. Spider is scheduled for flight over Aus- 
tralia in spring 2010 making measurements with 1600 detectors. The speaker em- 
phasised the important combined impact of these two experiments on determination 
of cosmological parameters in general. 

The ever increasing precision of CMB experiments warrants analysis of sub- 
tle observational effects inspired by the ideas from particle physics. The talk by 
Stephon Alexander "Can we probe leptogenesis from inflationary primordial bire- 
fringent gravitational waves" discussed a special mechanism of production of the 
lepton asymmetry with the help of cosmological birefringent gravitational waves. 
This mechanism was proposed in the recent paper-^^fil The mechanism assumes that 
gravitational waves are generated in the presence of a CP violating component of 
the inflaton field that couples to a gravitational Chern-Simons term (Chern-Simons 
gravity) . The lepton number arises via the gravitational anomaly in the lepton num- 
ber current. As pointed out by the speaker, the participating gravitational waves 
should lead to a unique parity violating cross correlation in the CMB. S. Alexander 
discussed the viability of detecting such a signal, and concluded by analyzing the 
corresponding observational constraints on the proposed mechanism of leptogenesis. 

Apart from being an arena for detecting relic gravitational waves, CMB is of 
course a powerful tool for cosmology in general. The final talk of the session, de- 
livered by Grant Mathews (co-authors D. Yamazaki, K. Ichiki and T. Kajino), 
discussed the "Evidence for a primordial magnetic field from the CMB tempera- 
ture and polarization power spectra" . This is an interesting subject as the magnetic 
fields are abundant in many astrophysical and cosmological phenomena. In partic- 
ular, primordial magnetic fields could manifest themselves in the temperature and 
polarization anisotropics of the CMB. The speaker reported on a new theoretical 
framework for calculating CMB anisotropics along with matter power spectrum in 
the presence of magnetic fields with power-law spectra. The preliminary evidence 
from the data on matter and CMB power spectra on small angular scales suggest 
an upper and a lower limit on the strength of the magnetic field and its spectral 
index. It was pointed out that this determination might be the first direct evidence 
of the presence of primordial magnetic field in the era of recombination. Finally, the 
author showed that the existence of such magnetic field can lead to an independent 
constraint on the neutrino mass. 

3. Analysis of WMAP data and outlook for Planck 

Along the lines of the presentation by Baskaran et al, we review the results of the re- 
cent analysis of WMAP 5-year data. The results suggest the evidence, although still 
preliminary, that relic gravitational waves are present, and in the amount character- 
ized by i? « 0.23, This conclusion follows from the likelihood analysis of WMAP5 
TT and TE data at lower multipoles £ < 100. It is only within this range of multi- 
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poles that the power in relic gravitational waves is comparable with that in density 
perturbations, and gravitational waves compete with density perturbations in gen- 
erating CMB temperature and polarization anisotropics. At larger multipoles, the 
dominant signal in CMB comes primarily from density perturbations. 

3.1. Likelihood analysis of WMAP data 

The analysis in paper^23124j based on proper specification of the likelihood func- 
tion. Since TT and TE data are the most informative in the WMAP5, the likelihood 
function was marginalized over the remaining data variables EE and BB. In what 
follows DJ^ and Dj^ denote the estimators (and actual data) of the TT and TE 
power spectra. The joint pdf for Df^ and Dj^ has the form 

f{Dr, DD = n^x^-^ {2i+".r^(|)(l - pl){<yjr{af)Y' 
J 1 / pgz X \ 1 

This pdf contains the variables (actual data) Df-^ {XX' — TT, TE) in quantities 
X = ti{DJ^ + NJ^) and z = nDj^ , where Nf^ is the temperature total noise 
power spectrum. The quantity n — {21 + l)/sky is the number of effective degrees of 
freedom at multipole where /sky is the sky-cut factor, /sky = 0.85 for WMAP and 
/sky = 0.65 for Planck. T is the Gamma-function. The quantities aj , af and pi, 
include theoretical power spectra Cf^^ and contain the perturbation parameters 
i?. As, ris and n*, whose numerical values the likelihood analysis seeks to find. 

The three free perturbation parameters i?. As, Us {rit = — 1) were determined 
by maximizing the likelihood function 

for £ = 2. ..100. The background cosmological model was fixed at the best fit ACDM 
cosmology.'23 The maximum likelihood (ML) values of the perturbation parameters 
(i.e. the best fit values), were found to be 

R = 0.229, Us = 1.086, As 1.920 x 10~^ (9) 

and nt = 0.086. The region of maximum likelihood was probed by 10,000 sample 
points using MCMC method. The projections of all 10,000 points on 2-dimensional 
planes R — Ug and R ^ Ag are shown in Fig. |3l 

The samples with relatively large values of the likelihood (red, yellow and green 
colors) are concentrated along the curve, which projects into approximately straight 
lines (at least, up to i? « 0.5): 

= 0.98 + 0.46i?, ^ = (2.27 - 1.53i?) x 10"^. (10) 

These combinations of parameters R, Us, As produce roughly equal responses in the 
CMB power spectra. The best fit model ([9]) is a particular point on these lines. 
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R = 0.229. The family of models ([TU)) is used in Sec l3.3l for setting the expectations 
for the Planck experiment. 

The marginalized 2-dimensional and 1-dimensional distributions are plotted in 
Fig. m and Fig. O respectively. The ML values for the 1-dimensional marginalized 
distributions and their 68.3% confidence intervals are given by 

R = 0.266 ± 0.171, Us = l.mt^^ml^ = x 10^^ (11) 



1.05 1,1 1,15 1.2 1.25 1.3 1.35 1 1.2 1.4 



Fig. 3. The projection of 10,000 samples of the 3-dimensional Kkehhood function onto the R — Us 
(left panel) and _R — As (right panel) planes. The color of an individual point in Fig. [3] signifies 
the value of the 3-dimensional likelihood of the corresponding sample. The black + indicates the 
parameters listed in Figure adopted from Zhao et al.I^ 



The derived results allow one to conclude that the maximum likelihood value 
of R persistently indicates a significant amount of relic gravitational waves, even 
if with a considerable uncertainty. The R = hypothesis (no gravitational waves) 
appears to be away from the R = 0.229 model at about la interval, or a little more, 
but not yet at a significantly larger distance. The spectral indices Ug, nt persistently 
point out to the 'blue' shape of the primordial spectra, i.e. > l,nt > 0, in the 
interval of wavelengths responsible for the analyzed multipoles £ < imax — 100. 
This puts in doubt the (conventional) scalar fields as a possible driver for the initial 
kick, because the scalar fields cannot support {3 > —2 and, consequently, cannot 
support Ug > l,7it > 0. 
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Fig. 4. The ML points (red x) and the 68.3% and 95.4% confidence contours (red solid lines) 
for 2-dimensionaI likelihoods: R — ris (left panel) and R — As (right panel). In the left panel, the 
WMAP confidence contours are also shown for comparison. Figure adopted from Zhao et 




1 1.2 1.4 1.6 1,8 2 2.2 2.4 2.6 



Fig. 5. 1-dimensional likelihoods for R (left), Us (middle) and As (right). Figure adopted from 
Zhao et al?^' 
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3.2. How relic gravitational waves can be overlooked in the 
likelihood analysis of data 

The results of this analysis differ from the conclusions of WMAP teamP^] xhe 



WMAP team found no evidence for gravitational waves and arrived at a 'red' spec- 
tral index Us = 0.96. The WMAP findings are symbolized by black dashed and blue 
dash-dotted contours in Fig. 21 It is important to discuss the likely reasons for these 
disagreements. 

Two main differences in the data analysis are as follows. First, the present anal- 
ysis is restricted only to multipoles £ < 100 (i.e. to the interval in which there 
is any chance of finding gravitational waves), whereas the WMAP team uses the 
data from all multipoles up to i ^ 1000, keeping spectral indices constant in the 
entire interval of participating wavelengths (thus making the uncertainties smaller 
by increasing the number of included data points which have nothing to do with 
gravitational waves). Second, in the current work, the relation nt — ria — 1 implied 
by the theory of quantum-mechanical generation of cosmological perturbations is 
used, whereas WMAP team uses the inflationary 'consistency relation' r — ~8nt, 
which automatically sends r to zero when nt approaches zero. 

It is important to realize that the inflationary 'consistency relation' 



is a direct consequence of the single contribution of inflationary theory to the subject 
of cosmological perturbations, which is the incorrect formula for the power spectrum 
of density perturbations, containing the "zero in the denominator" factor: 



pendently of the curvature of space-time and strength of the generating gravitational 
field characterized by H. To make the wrong theory look "consistent" , inflationary 
model builders push H/Hpi down, whenever e goes to zero (for example, down to 
the level marked by the minutely inclined line "legitimate transition in inflation- 
ary theory" in Fig. [T]), thus keeping at the observationally required level and 
making the amount of relic gravitational waves arbitrarily small. In fact, the most 
advanced inflationary theories based on strings, branes, loops, tori, monodromies, 
etc. predict the ridiculously small amounts of gravitational waves, something at the 
level of r « 10^^^, or so. [There is no doubt, there will be new inflationary theories 
which will predict something drastically different. Inflationists are good at devising 
theories that only mother can love, but not at answering simple physical questions 
such as quantization of a cosmological oscillator with variable frequency. For a more 
detailed criticism of inflationary theory, see Grishchuk.-j Obviously, the analysis in 
paper^^^'^il ^joes not use the inflationary theory. 



r — 16e — —Sut 




The "zero in the denominator" factor e is e = —H/H^. This factor tends to zero in 
the limit of standard de Sitter inflation H — (any horizontal line in Fig. [T]) inde- 
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Baskaran et al concluded that in the conditions of relatively noisy WMAP data 
it was the assumed constancy of spectral indices in a broad spectrum^'' that was 
mostly responsible for the strong dissimilarity of data analysis results with regard 
to gravitational waves. Having repeated the same analysis of data in the adjacent 
interval of multipoles 101 < £ < 220, the authors of papei^^^ found a somewhat 
different (smaller) spectral index n,, in this interval. They came to the conclusion 
that the hypothesis of strictly constant spectral indices (perfectly straight lines in 
Fig. [1]) should be abandoned. It is necessary to mention that the constancy of Us 
over the vast region of wavenumbers, or possibly a simple running of Ug, is a usual 
assumption in a number of works. 1^31221 

It is now clear why it is dangerous, in the search for relic gravitational waves, 
to include data from higher multipoles, and especially assuming a strictly constant 
spectral index Ug. Although the higher multipole data have nothing to do with 
gravitational waves, they bring Ug down. If one postulates that Ug is one and the 
same constant at all £'s, this additional 'external' information about Ug affects 
uncertainty about R and brings R down. This is clearly seen, for example, in the 
left panel of Fig. E) The localization of Ug near, say, the line ng = 0.96 would cross 
the solid red contours along that line and would enforce lower, or zero, values of 
R. However, as was shown,'^^! Ug is sufficiently different even at the span of the 
two neighbouring intervals of ^'s, namely 2 < £ < 100 and 101 < i < 220. These 
considerations, as for how relic gravitational waves can be overlooked in a data 
analysis, have general significance and will be applicable to any CMB data. 

3.3. Forecasts for the Planck mission 

The final part of the presentation by Baskaran et al dealt with forecasts for the 
recently launched Planck satellite.^ The ability of a CMB experiment to detect 
gravitational waves is characterized by the signal-to-noise ratio defined bj^^SMl 

S/N = R/AR, (12) 

where the numerator is the "true" value of the parameter R (its ML value or the 
input value in a numerical simulation) while AR in the denominator is the expected 
uncertainty in determination of R from the data. 

In formulating the observational expectations, all of the available information 
channels (i.e. TT, TE, EE and BB correlation functions) were taken into account. 
The calculations were performed using the Fisher matrix formalism. The total un- 
certainty entering the Fisher matrix calculations is comprised of instrumental and 
foreground noises^ 24i[30] ^^^1 as the statistical uncertainty of the inherently ran- 
dom signal under discussion. The possibility of partial removal of contamination 
arising from foreground sources was quantified by the residual noise factor o-fs. The 
three considered cases included a^^ = 1 (no foreground removal), a^^ = 0.1 (10% 
residual foreground noise) and a^^ = 0.01 (1% residual foreground noise). In or- 
der to gauge the worst case scenario, the 'pessimistic' case was added. It assumes 
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a^^ = 1 and the instrumental noise at each frequency Vi four times larger than the 
nominal value reported by the Planck team. This increased noise is meant to mimic 
the situation where it is not possible to get rid of various systematic effectsj^ the 
E-B mixingf22l cosmic lensingj^ etc. which all affect the BB channel. 

The total S/N for one parameter family of models defined by the large WMAP5 
likelihoods pH)) is plotted in Fig. [51 This graph is based on the conservative as- 
sumption that all unknown parameters R, Hs and Ag are evaluated from one and 
the same dataset. Four options are depicted: cr^^ — 0.01,0.1,1 and the pessimistic 
case. The results for the ML model ([9|) are given by the intersection points along the 
vertical line R = 0.229. Specifically, S/N = 6.69, 6.20, 5.15 for cr*s = 0.01, 0.1, 1, 
respectively. The good news is that even in the pessimistic case one gets S/N > 2 for 
R > 0.11, and the Planck satellite will be capable of seeing the ML signal R ~ 0.229 
at the level S/N = 3.65. 




It is important to treat separately the contributions to the total signal-to-noise 
ratio supplied by different information channels and different individual multipoles. 
The (S/N)^ calculated for three combinations of channels, TT + TE + EE + BB, 
TT + TE + EE and BB alone, is shown in Fig. [71 Four panels describe four noise 
models: a^^ — 0.01, 0.1, 1 and the pessimistic case. Surely, the combination TT + 
TE + EE + BB is more sensitive than any of the other two, i.e. TT + TE + EE and 
BB alone. For example, in the case a^^ =0.1 the use of all correlation functions 
ensures S/N which is ~ 50% greater than BB alone and 30% greater than 
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TT + TE + EE. The situation is even more dramatic in the pessimistic case. The 
ML model © can be barely seen through the S-modes alone, because the BB 
channel gives only S/N = 1.75, whereas the use of all of the correlation functions 
can provide a confident detection with S/N = 6.48. Comparing TT + TE + EE 
with BB one can see that the first method is better than the second, except in the 
case when a^^ = 0.01 and R is small {R < 0.16). In the pessimistic case, the role of 
the BB channel is so small that the TT + TE + EE method provides essentially 
the same sensitivity as all channels TT + TE + EE + BB together. 

Finally, Fig. [8] shows the contributions to the total signal-to-noise ratio from 
individual multipoles. It can be seen that a very deep foreground cleaning, a^^ = 

0. 01. makes the very low (reionization) multipoles £ 10 the major contributors to 
the total (S/N)'^, and mostly in the BB channel. However, for large a^^ — 0.1, 1, 
and especially in the pessimistic case (see the lower right panels in Fig. [8]), the 
ability of the BB channel becomes very degraded at all ^'s. At the same time, as 
Fig. E] illustrates, the ^-decomposition of (S/N)^ for TT + TE + EE combination 
depends only weakly on the level of a^^. Furthermore, in this method, the signal 
to noise curves generally peak a.t £ ^ (20 — 60). Therefore, it will be particularly 
important for Planck mission to avoid excessive noises in this region of multipoles. 

4. Conclusions 

In general, the OCl parallel session was balanced and covered both theoretical and 
experimental aspects of primordial gravitational waves and the CMB. Thanks to 
excellent contributions of all participants, this subject received a new momentum. 
Hopefully, new observations and theoretical work will bring conclusive results in the 
near future. 
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